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Abstract 

Stars twinkle to the eye through atmospheric turbulence, but planets, because of 
their larger angular size, do not. Similarly, scintillation due to the local interstel¬ 
lar medium will modulate the radio flux of gamma-ray-burst afterglows and may 
permit indirect measurements of their angular sizes. The amplitude of refractive 
scintillation is of order ten percent at ten gigahertz unless the source size is much 
larger than the expected size, of order ten microarcseconds. Diffractive scintillation 
is marginally possible, depending sensitively on the source size, observing frequency, 
and scattering measure of the interstellar medium. 
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1 INTRODUCTION 


At the time of writing, three gamma-ray bursts discovered by the BeppoSAX 
satellite have been associated with transient sources at longer wavelengths 
([Costa et al. 1997], [Feroci et al. 1997], [Piro et ah 1997], [Costa et ah 1997], 
[Bond 1997]). Absorption lines of iron and magnesium at z = 0.835 in the op¬ 
tical counterpart of GRB0508 appear to prove that this burst, at least, lies at 
a cosmological distance ([Metzger et ah 1997]). GRB0508 is also the first burst 
associated with a radio transient ([Frail & Kulkarni 1997], [Pooley & Green 1997]). 


These Xray-to-radio “afterglows” are broadly consistent with models based 
on relativistic blastwaves at cosmological distances ([Rees & Meszaros 1992], 
[Paczynski & Rhoads 1993], [Waxman 1997], [Wijers et ah 1997], and refer¬ 
ences therein). In these models, the shock Lorentz factor is ^ 300 at the 
time of the gamma-ray event and falls thereafter as r s oc t -3 / 8 , t being the 
time measured at Earth. The bulk Lorentz factor of the postshock mate¬ 
rial is r ps = 2 - 1 / 2 r s . For synchrotron emission, the flux density F v at fre¬ 
quency v is predicted to peak at a time t v oc z/~ 3//2 if the source is opti¬ 
cally thick ([Paczynski & Rhoads 1993]), or t v oc z/~ 2 / 3 if optically thin (cf. 
[Waxman 1997], [Wijers et ah 1997]). 


A lower limit to the angular size 9 S = r s /D of the radio source follows if 
the brightness temperature is not larger than the Compton limit for inco¬ 
herent synchroton radiation, Tb <; r ps (l + /3 ps )T max with T max « 10 12 K 
([Kellermann & Pauliny-Toth 1969]): 


0 S > 1.5 


1 + z 

1/2 ( F„ 

_(i + /3 P s)r ps _ 

V mJ y 


1/2 


V 


— 1/2 ( -Rnax \ 

10 V10 12 A / 


- 1/2 


/i cl.S. ? 


( 1 ) 


where iq 0 = z^/(10 GHz). [Sari (1997)] has pointed out that the radius of the 
shock at observed time t is R s = 8r 2 ct for T s 1 since T s oc t -3 / 8 . Most of 
the emission is seen from a disk of radius R s /T ps , so 


= Ct f 8V2r s r s > 1; 

d \ (5/2)/3 s (3 S = (1 - T~ 2 ) 1 / 2 < 1 


( 2 ) 


where d , the proper motion distance of the source, is (1 + ^) times the angular- 
size distance (cf. [Weinberg 1972]), and the nonrelativistic formula above fol¬ 
lows from the Sedov solution. Comparison of relations (1) and (2) yields 
lower bounds on T s as well as 0 S . Thus for example, the radio counterpart 
of GRB0508 had a flux density of 0.61 ± 0.04 mJy at 8.46 GHz 6.2 days af¬ 
ter the burst ([Frail & Kulkarni 1997]). Using a smooth interpolation between 
the ultrarelativistic and nonrelativistic formulae above, and for the relation 
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between T s and r ps , we find / 3 S ^ 0.5 and 9 S ^ 1.5 p.s. at the time of this ob¬ 
servation, assuming z = 1 and d = 10 28 cm (3 Gpc). The spherical blastwave 
models cited above predict somewhat larger values, depending on parameters. 
For example, [Sari (1997)]’s estimates imply T s ~ 2.6(E 52 /rii) 1 ^ s at t = 6.2 d, 
where E$ 2 is the energy of the blastwave in units of 10 52 erg and ri\ is the den¬ 
sity of the external medium in cm -3 ; this implies 9 S « 9 ju a.s. for the same 
values of d and z as above. “Jet” models, in which the ejecta are t unn eled into 
a small solid angle, predict angular sizes ~ 4 times smaller than corresponding 
spherical ones but may already be ruled out because they vastly underpredict 
the observed radio flux ([Rhoads 1997]). 

Thus, if GRB0508 proves to be typical, it is unlikely that radio afterglows will 
be directly resolved by VLB1. An upper limit of 300 fi a.s. has been set by the 
VLBA eight days after the burst ([Taylor, Beasley, & Frail, 1997]). 

Afterglows may, however, be resolved indirectly by interstellar scattering. In¬ 
homogeneities in the density of free electrons (N e ) affect the refractive index 
(n) governing the propagation of radio waves through the interstellar medium: 
5n = — (r e \ 2 /2ir)5N e , where r e = e 2 /m e c 2 , and A is the wavelength. Exten¬ 
sive evidence supports the conclusion that these inhomogeneities span a broad 
range of spatial scales (at least 10 9 to 10 14 cm) with a power-law power spec¬ 
trum 

®nM = ( 27r )' 3 / d 3 s(6N e (x)5N e (x + s))e^ = C 2 N q~ u/3 , (3) 


where (f is the spatial wavenumber conjugate to separation s and is measured 
in radians per unit length (cf. [Armstrong, Rickett, & Spangler 1995]). Small- 
scale atmospheric density fluctuations have power spectra of the same math¬ 
ematical form, so that radio scintillation has much in common with optical 
seeing ([Coulman 1985]). 

Scintillation modulates the flux in two conceptually different ways (cf. [Rickett 1990]): 

(i) Refractive scintillation, the random focusing and defocusing of rays, can be 
understood and analyzed entirely within geometric optics. Given a turbulent 
spectrum such as (3), the strength of refractive scintillation varies smoothly 
with the angular size of the source (inhomogeneities smaller than the projected 
source size have no effect because surface brightness is conserved) and with 
wavelength (5n oc A 2 ). Refractive scintillation is insensitive to the observing 
bandwidth. 

(ii) Diffractive scinti 1 lation is a physical-optics effect; it can be explained as the 
interference among multiple paths from source to receiver. Since it requires the 
interfering rays to be mutually coherent, diffractive scintillation is quenched 
when the angular size of the source or the bandwidth of observation is too 
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large. Diffractive scintillation leads to a Rayleigh distribution of flux, with 
fluctuations equal to the mean: P(F) = exp(— F/(F)), (F 2 ) = 2(F) 2 . 

Among common radio sources, only pulsars are compact enough to display 
diffractive scintillation, but refractive scintillation affects interstellar masers 
and extragalactic sources as well. 

In §2, we consider diffractive scintillation of radio afterglows. We find that 
this is just possible but very sensitive to the actual angular size of the af¬ 
terglow and to the radio frequency of observation. Detection of diffractive 
scintillation would strongly constrain the Lorentz factor of the burst. In §3, 
we describe the more robust predictions for refractive scintillation, which can 
also be used to constrain the angular size. Flux variations relative to the mean 
~ 0.1(# s / /i a.s.) -7 ' /6 z/ ] j ) 2 are predicted on timescales ~ 10 hr. A summary and 
some suggestions for observing strategies are presented in §4. 


2 DIFFRACTIVE SCINTILLATION 


For convenience, we imagine in this section only that the scattering medium 
is compressed onto a plane perpendicular to the line of sight at distance z sc 
(“thin-screen approximation”). The characteristic deflection due to scattering 

is 
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Here SM = / C%(z)dz, where z is a coordinate along the line of sight, is called 
the scattering measure. We have scaled it by a value typical for extragalactic 
sources at high galactic latitudes ([Spangler et al. 1993]). 


The precise definition of 9 d involves the phase structure function, D^s), which 
is the mean-square phase difference accumulated along parallel lines of sight 
separated by transverse distance s. It follows from eq. (3) that D^s) = 
( S / Sd ) 5 / 3 for an appropriately chosen constant s d . Then 9 d = (/cs d ) _1 , where 
k = 27t/A is the wavenumber. Apart from a numerical factor, s d is equivalent 
to the Fried parameter r d of optical seeing. 


The scatter-broadened image acts somewhat as a Michelson stellar interfer¬ 
ometer, with a characteristic spacing between the arms a ~ z sc 9 d . Fringes 
develop on the observer’s plane if the angular size of the source is smaller 
than the resolution of the interferometer, 9 a = (/ca) -1 . Diffractive scintillation 
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is caused by the motion of the observer through this fringe pattern. Therefore, 
diffractive scintillation is seen only if 

9 S < (kz sc $ d ) —1 = 2.25z/® 0 /5 ^ kpc (S , M_3. 5 )' 3/5 /za.s., (5) 


where SM_ 3 5 = SM/{ 10 3 ' 5 m 20//3 kpc). Condition (5) can be satisfied if the 
source size is close to the lower limit discussed in §1. 

In order to see diffractive flux variations of relative order unity, however, 
it is also necessary that D^r-p) > 1, or equivalently that Sd < ry, where 
ry = yjd scr /k is the Fresnel length (cf. [Rickett 1990]). The rays of geometric 
optics have an intrinsic width ~ ry, and in order to create interference effects, 
the phase fluctuations must be strong enough to create distinct rays. If this 
condition is not fulfilled, there is only weak scintillation (weak modulation of 
the flux): a good optical site is in the weak-scintillation regime, so that the 
stars overhead do not twinkle. The condition for strong scintillation in the 
present case is 

v < 10.4(SM_3. 5 ) 6/17 i / , 1 V GHz - (6) 

At the maximum frequency allowed by eq. (6), the source-size limit becomes 

9 s (u max ) < 2 . 35 ( SM _ 3. 5 )- 3/1 Vci'kp" Ra.s.. (7) 

Intriguingly, this is comparable to the minimum source size estimated in §1 
on physical grounds. Both the scattering measure SM and the effective dis¬ 
tance z scr scale as esc b with galactic latitude [6(GRB0508) ~ 27°]. Hence the 
maximum measurable size (7) is oc (sinfr) 14 / 17 . 

Future afterglows may therefore be indirectly resolvable by diffractive scintil¬ 
lation when they occur at moderately high latitude, especially if they can be 
detected shortly after the burst 1 d). Scintillation will be quenched at some 
frequency below (6) because of the z/ 6 / 5 dependence of the constraint (5), and 
the transition frequency will provide an indirect measure of the source size. 
The scattering measure along the line of sight can be internally estimated from 
the scintillation timescale, 

tdm = — = 3.1zqo 5 (SM_ 3 . 5 )“ 3/5 ( — U± A hr. (8) 

v± V30 km s / 

In very strong scintillation, diffractive flux variations are uncorrelated at radio 
frequencies differing by more than the decorrelation bandwidth 

A^c = » 7.64 2/ AM-3,r«-A-I. kpc GHz. (9) 

d. scr 
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As this formula shows, the frequencies of interest are not far from the bor¬ 
derline between strong and weak scintillation, where Az/<j c ~ v- On the one 
hand, this is fortunate because the afterglows are faint and must be observed 
with a broad bandwidth. On the other hand, it is unfortunate because 
measures an independent combination of SM and d screen . 

The angular size limit (5) and the decorrelation bandwidth (9) become more 
severe as the distance to the scattering screen increases, and it might therefore 
be supposed that scattering in the gamma-ray burst’s host galaxy (if it exists) 
or in the intergalactic medium (IGM) would suppress diffractive scintillation. 
At worst, however, such scattering has the effect of replacing one incoherent 
source (of size 9 S ) with an effective incoherent source of somewhat larger an¬ 
gular size 9 S . Diffractive scintillation clue to the Galactic interstellar medium 
(ISM) will still occur provided that 9 S satisfies the limit (5) with d scr deter¬ 
mined by the path length through the Galactic ISM. Also, a deflection angle 
Od in the host galaxy contributes to the effective source size in proportion to 

AMe d ^ io-«9 d . 

dhost 


Hence, scattering in the host can probably be neglected, for the same reason 
that orbiting satellites looking down enjoy much better seeing than ground- 
based telescopes looking up. 

Scattering by the IGM is negligible presuming (in the absence of better infor¬ 
mation) that the scattering measure scales as the path length times the square 
of the mean electron density: 

SMiqm 1 Gpc /2 x 1(T 7 cm _3 \ 2 4 

AMism 1 kpc \ 0.02 cm -3 ) 


We have assumed a fully-ionized IGM containing most of the baryons allowed 
by primordial nucleosynthesis (fibh| 0 = 0.05±0.01: [Izotov, Thuan, & Lipovetsky 1997]) 
and taken the local electron density from [Taylor & Cordes 1993]. In the case 
of GRB0508, scattering by an intervening galaxy associated with the z = 

0.835 absorption system might be comparable to that due to the local ISM, 
but there appears to be no intrinsically bright galaxy on the line of sight 
[[Fruchter & Bergeron 1997]], so it is likely that the absorption comes from 
the outer halo of a galaxy, where the electron density and scattering measure 
are likely to be small. 
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3 REFRACTIVE SCINTILLATION 


Following [Coles et al. 1987] (henceforth CFRC), the normalized spatial cor¬ 
relation of the flux due to refractive scintillation is 
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Here the line of sight runs from the observer at z = 0 to the source at (ef¬ 
fectively) z — oo. The electron-density spectrum <3>jv e [eq. (3)] depends on z 
via the coefficient C% —> Cjf(z), which reflects changes in the strength of the 
turbulence along the line of sight. We take C%{z) = C%{$) exp[— (zsinb/H) 2 ], 
where Cfj(0) is the local value near the Sun, and z sin b is the height above the 
plane. The flux correlation depends upon the average (or sum) of the phase 
fluctuations along the line of sight, hence the z component of q is set to zero 
in eq. (10). 


The quantity V(r) is the visibility of the source on baseline r, which like q 
and s in this formula, is a vector transverse to the line of sight. We assume a 
gaussian image brightness distribution for the source, /„(#) oc exp(— 9 2 /29 2 ), 
so that V(qz/k) = exp(—q 2 z 2 9 2 /2). 


Finally, 


D'Jr, z) = A7rr 2 e X 2 J d 2 q$ Ne (q x , q y \ q z = 0 , z) 

= 2 4 /3 3?r2 r(7/6) 2^ ( } 2 2 5/3 
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is the differential phase structure function. In other words, D'Jr, z)dz is the 
contribution of the slab (z, z+dz ) to the mean-square phase difference between 
lines of sight separated by r. The quantity z< = min( 2 , z). 


The three hnal factors in equation (10) can be regarded as the squares of 
“visibilities” due to the intrinsic source size, to the scattering, and to Fresnel 
optics. The characteristic angular scales associated with these visibilities at 
z = H esc b are 9 S , 9 d , and the Fresnel angle 

9y = {kHescb)~ 1 ^ 2 — 2.fia.s.. ( 12 ) 

V k P c J 
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Actually, the Fresnel visibility at z — Hcscb is 


V F (r) 


sm(k6 F r/2) 

k0 F r/2 


( 13 ) 


so that Vp(0) = 1. All three visibilities have the effect of restricting the range of 
wavenumbers q that contribute to C(s). They occur multiplicatively in equa¬ 
tion (10) because the effective source size seen by the observer is a convolution 
of the intrinsic surface brightness distribution I v {6) with the scatter-broadened 
image of a point source, and with a sort of Fresnel image. 

The root-mean-square variation in the flux relative to its mean value is \JC{ 0) = 
m R , the modulation index. From equation (10), 
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This extends the results of CFRC, who gave explicit formulae for m R in the 
limit of large intrinsic size, 6 S 0 dl 0 F . Our result (14) is an interpolation 
formula. The numerical coefficients in the definition of 0 cF have been chosen 
so that equation (14) is exact when any one of the three characteristic angles 
(# s , 0 dl 0 F ) is much larger than the other two. 

Refractive scintillation occurs on the timescale required for the line of sight 
to cross the image at a mean distance d scr ~ H esc b/2, 

t,„ = e ‘ sHaXb = 7.0 ( ( V1 y hr 

2v± y 10/ia.s./ \ kpc / V30 km s 1 / 

assuming that the motion of the line of sight is more rapid than the internal 
velocities of the medium. The full temporal correlation function is then C(v±r) 
at lag r. Note that unlike t di g [eq. (8)], t re f remains constant or decreases with 
increasing frequency, depending on the intrinsic image size. At frequencies 
above limit (6), where the scintillation is weak, the two timescales merge into 
a single one: the time required for the line of sight to cross the Fresnel length. 

In principle, it is possible to estimate the angular size from the amplitude 
and timescale of refractive scintillation at a single frequency. In practice, it 




may be difficult to disentangle refractive flux variations from instrumental 
and atmospheric noise; futhermore, the parameters C'jy-(O) and H are imper¬ 
fectly known. Therefore it will be safer to measure mu at several frequencies. 
Equations (14), (15), and (4) imply that m r oc z / -2 at relatively high frequen¬ 
cies, where 0 e s ~ 9 S is frequency-independent; and that rriR oc v +l ‘/ m a t low 
frequencies, where 0 e s ~ 6& oc z/ -11 / 5 . Hence 771 r(z/) has a peak: 
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These formulae assume that 0 S is substantially larger than the Fresnel angle 
(12) at the peak. 


4 Discussion 


If radio afterglows of gamma-ray bursts are indeed produced by relativistic 
blastwaves at cosmological distances, then their small angular size and high 
surface brightness puts them in an extremely interesting part of parameter 
space with respect to interstellar scintillation. Both diffractive and refractive 
scintillation are possible for these sources. 

Diffractive scintillation is sensitive to the radio frequency of observation, u, 
and the angular size of the source, 6 S . It may have been marginally possible 
to have seen diffractive scintillation in GRB0508 at v & 10 GHz if 0 S £ 
2 fi a.s. [eq. (5)]. Future afterglows from sources at high galactic latitude will 
present other opportunities for observing diffractive scintillation, which will 
be recognized because 

• the flux variations are of unit strength, (F 2 ) = 2(F) 2 [but contrary to 
the situation with pulsars, scintillation will be correlated accross a broad 
frequency band, cf. eq. (9)]; 

• the timescale decreases with decreasing radio frequency, and is on the order 
of a few hours [eq. (8)]; 

• the effect vanishes abruptly below a frequency determined by the scattering 
measure and the angular size of the source [eq. (6)]. 

Diffractive scintillation will allow the interstellar medium to be used as a 
Michelson stellar interferometer to resolve angular sizes measured in micro- 
arcseconds, far too small for conventional VLBI. 
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Refractive scintillation will occur at all frequencies and for all angular sizes, 
though its amplitude depeneds algebraically on both of thes parameters. The 
timescale will range upward from a few hours. The amplitude, though less 
than unity, will be large (~ 10%) unless 9 S , and hence the Lorentz factor of 
the blastwave, are much larger than predicted by current models [eq. (14)]. Re¬ 
fractive scintillation therefore can also be used to constrain 9 S and T s , though 
it will not provide quite so sharp a test as diffractive scintillation. 

Some implications for observing strategies are the following: 

One wants to collect data at several frequencies, roughly in the range 2 — 
20 GHz. 

For diffractive scintillation, the flux should be sampled on timescales ~ 1 hr, 
and the first data points should be taken as early as the sensitivity of the 
equipment permits, since the source will be smaller and less likely to be over¬ 
resolved by the scattering [see eqs. (5) & (7)]. One should probably attempt 
to see diffractive scintillation in sources at high galactic latitude 30°) only, 
once again because these are least likely to be overresolved. Ideally one would 
wait for a source close to a lines of sight on which the scattering parame¬ 
ters have been well determined from observations of extragalaetic sources and 
globular-cluster pulsars. 

For refractive scintillations, many of the same considerations apply. In this 
case, however, one can afford to make observations later in the lifecycle of the 
afterglow, at longer time intervals, and at lower galactic latitude. It will be 
important, however, to minimize or accurately subtract the contribution of 
receiver and atmospheric noise to the flux variations, since the amplitude of 
the modulation index is critical to the estimate of the angular size. 

I am indebted to Michael Rupen, Andrew Ulmer, Michael Richmond, Bohdan 
Paczynski, and the rest of the Observatory coffee klatsch for stimulating dis¬ 
cussions and for criticism of the first draft. Special thanks are owed to Michael 
Richmond for help with the IAU circulars and for his ephemeris software. 
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